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Abstract. Consider the random normal matrix ensemble associated with a potential on the plane 
which is sufficiently strong near infinity. It is known that, to a first approximation, the eigenval- 
ues obey a certain equilibrium distribution, given by Frostman's solution to the minimum energy 
problem of weighted logarithmic potential theory. On a finer scale, one can consider fluctuations of 
eigenvalues about the equilibrium. In the present paper, we give the correction to the expectation of 
fluctuations, and we prove that the potential field of the corrected fluctuations converge on smooth 
test functions to a Gaussian free field with free boundary conditions on the droplet associated with 
the potential. 



Given a suitable real "weight function" in the plane, it is well-known how to associate a 
corresponding (weighted) random normal matrix ensemble (in short: RNM-ensemble). Under 
reasonable conditions on the weight function, the eigenvalues of matrices picked randomly from 
the ensemble will condensate on a certain compact subset S of the complex plane, as the order of 
the matrices tends to infinity. The set S is known as the droplet corresponding to the ensemble. 
It is well-known that the droplet can be described using weighted logarithmic potential theory 
and, in its turn, the droplet determines the classical equilibrium distribution of the eigenvalues 
(Frostman's equilibrium measure). 

In this paper we prove a formula for the expectation of fluctuations about the equilibrium 
distribution, for linear statistics of the eigenvalues of random normal matrices. We also prove 
the convergence of the potential fields corresponding to corrected fluctuations to a Gaussian free 
field on S with free boundary conditions. 

Our approach uses Ward identities, that is, identities satisfied by the joint intensities of the 
point-process of eigenvalues, which follow from the reparametrization invariance of the partition 
function of the ensemble. Ward identities are well known in field theories. Analogous results in 
random Hermitian matrix theory are known due to Johansson |13], in the case of a polynomial 
weight. 



General notation. By D(a, r) we mean the open Euclidean disk with center a and radius r. By 
"dist" we mean the Euclidean distance in the plane. If A„ and B„ are expressions depending on 
a positive integer n, we write A„ < B n to indicate that A n < CB„ for all n large enough where 
C is independent of n. The notation A n x £>„ means that A n < £>„ and B n < A n . When p. is a 
measure and / a ,u-measurable function, we write ju(f) = J f d\i. We write d = j(d/dx - id/dy) 
and d = j(d/dx + id/dy) for the complex derivatives. 
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1. Random NORMAL MATRIX ENSEMBLES 

1.1. The distribution of eigenvalues. Let Q : C — » RU {+00} be a suitable lower semi-continuous 
function subject to the growth condition 

(1.1) liminf-^- > 1. 

z^oo log |z| 

We refer to Q as the weight function or the potential. 

Let N„ be the set of all n x n normal matrices M, i.e., MM* = M*M. The partition function on Af„ 
associated with Q is the function 



where dM„ is the Riemannian volume form on N„ inherited from the space C" of all n X n matrices, 
and where trace Q : N n —* K U is the random variable 

trace Q(M)= £ Q(Ay), 

Aj e spec(M) 

i.e., the usual trace of the matrix Q(M). We equip N„ with the probability measure 

dPn = -L e -2» trace Q(M) ^ 

and speak of the random normal matrix ensemble or "RNM-ensemble" associated with Q. 

The measure P„ induces a measure P„ on the space C" of eigenvalues, which is known as the 
density of states in the external field Q; it is given by 



Here we have put 



dP„(A) = -L e - H " (A) dA n (A), A = (Aj)1 e C". 



H„(A) = J" log — ^— 1 +2nf Q(Aj), 
^ k |Ay - A fc | ^ 



and dA„(A) = d 2 Ai • • • rf 2 A„ denotes Lebesgue measure in C", while Z„ is the normalizing constant 
giving P„ unit mass. By a slight abuse of language, we will refer to Z n as the partition function of 
the ensemble. 

Notice that H n is the energy (Hamiltonian) of a system of n identical point charges in the plane 
located at the points Ay, under influence of the external field 2nQ. In this interpretation, P„ is 
the law of the Coulomb gas in the external magnetic field 2nQ (at inverse temperature /3 = 2). In 
particular, this explains the repelling nature of the eigenvalues of random normal matrices; they 
tend to be very spread out in the vicinity of the droplet, just like point charges would. 

Consider the n-point configuration ("set" with possible repeated elements) {Ay}" of eigenvalues 
of a normal matrix picked randomly with respect to f n . In an obvious manner, the measure F„ 
induces a probability law on the n-point configuration space; this is the law of the n-point process 
= {Ay}™ associated to Q. 

It is well-known that the process W n is determinantal. This means that there exists a Hermitian 
function K„, called the correlation kernel of the process such that the density of states can be 
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represented in the form 

dP n (A) = i det (k„(A ; , hjf dAM), AeC". 

One has 

K n {z,w) = K n (z,w)e- n{Q{z)+Q{w)) / 
where K n is the reproducing kernel of the space P„ (e~ 2n< ^ of analytic polynomials of degree at 
most n — 1 with norm induced from the usual L 2 space on C associated with the weight function 
g -2nQ Alternatively, we can regard K„ as the reproducing kernel for the subspace 

W n = \pe~ n ®; p is an analytic polynomial of degree less than n) c L 2 (C). 
We have the frequently useful identities 

f(z)= f f(w)Kn(z,w) tfw, feW n , 
Jc 

and 

K„(z,z) d 2 z = n. 

Moreover, if E„ denotes expectation with respect to P,„ and if / is a function depending on k 
variables, where k < n, then 

En [/(Ax, . . . , A k )] = (n ~ n f M B k [/(Aj A k )] . 

We refer to 0, (TBI, JU, lHOl , O, [14] for more details on point-processes and random matrices. 

1.2. The equilibrium measure and the droplet. We are interested in the asymptotic distribution 
of eigenvalues as n, the size of the matrices, increases indefinitely. Let u n denote the one-point 
function of P,„ i.e., 

«n(A) = -K„(A,A), AeC. 
n 

With a suitable function / on C, we associate the random variable Tr„ [/] on the probability space 
(C",P„)via 

n 

Tr„[/](A) = £/(A,). 

i=i 

The expectation is given by 

E n (Tr„[/]) = n I /■«„. 
Jc 

According to Johansson (see 1 10]) we have weak-star convergence of the measures 

da„(z) = u n (z)d 2 z 

to some probability measure a = a(Q) on C. 

In fact, a is the Frostman equilibrium measure of the logarithmic potential theory with external 
field Q. We briefly recall the definition and some basic properties of this probability measure, cf . 
1 16 1 and |10| for proofs and further details. 

Let S = supp a and assume that Q is C 2 -smooth in some neighbourhood of S. Then S is 
compact, Q is subharmonic on S, and a is absolutely continuous with density 

1 .~ 
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We refer to the compact set S = Sq as the droplet corresponding to the external field Q. 

Our present goal is to describe the fluctuations of the density field \i n = L/=i ^a, around the 
equilibrium. More precisely, we will study the distribution (linear statistic) 

/ » ti n (J) - no(f) = Tr„[/] - na{f), f e C °°(C). 

We will denote by v„ the measure with density n(u n - u), i.e., 

v n [f] = E„ [Tr„[/]] - na{f) = n{a n - a)(f), f 6 C~(C). 

1 .3. Assumptions on the potential. To state the main results of the paper we make the following 
three assumptions: 

(Al) (smoothness) Q is real analytic (written Q e C) in some neighborhood of the droplet S = Sq; 

(A2) (regularity) AQi=0in S; 

(A3) (topology) dS is a C'-smooth Jordan curve. 

We will comment on the nature and consequences of these assumptions later. Let us denote 

L = log AQ. 

This function is well-defined and C a ' in a neighborhood of the droplet. 

1.4. The Neumann jump operator. We will use the following general system of notation. If g is 
a continuous function defined in a neighborhood of S, then we write g s for the function on the 
Riemann sphere C such that g s equals g in S while g s equals the harmonic extension of g \ gs to 
C \ S on that set. 

If g is smooth on S, then 

Nag := Q := int(S), 

where n is the (exterior) unit normal of Q. We define the normal derivative A/q»£ for the 
complementary domain Q, := C \ S similarly. If both normal derivatives exist, then we define 
(Neumann's jump) 

Ng = N ds := Nag + N a ,g. 
By Green's formula we have the identity (of measures) 
(1.2) Ag s =Ag-l a + N(g s )ds, 

where ds is the arclength measure on dS. 

We now verify l|1.2t . Let <p be a test function. The left hand side in (L2) applied to <p is 

f <pAg s = f g s Acp= fgAcp+ f g s Acp, 
Jc Jc Js Jc\s 

and the right hand side is 

f<t>Ag + j ct>N(g s )ds. 

Thus we need to check that 

f (gA<p - QAg) + f (g s Acp - tpAg 5 ) = f ct>N(g s )ds. 
Js Jc\s J 



RANDOM NORMAL MATRICES AND WARD IDENTITIES 



5 



But the expression in the left hand side is 

f(gd n (p - (f>d„g)ds + f(gd m (f> - (pd m g s )ds = - J((pd n g + (pd m g s )ds = J (f> N (g s )ds, 
and l|1.2t is proved. 

1.5. Main results. We have the following results. 
Theorem 1.1. For all test functions f e C^°(C), the limit 

v(f) := lim v„(f) 

n—>oa 

exists, and 



f A/ + T/AL+ f fN(L s )ds 

JS JS JdS 



Equivalently, we have 



v«-v=i-A(l s+ L s ) 
in the sense of distributions. 



Theorem 1.2. Let h 6 C^(C) be a real-valued test function. Then, asn — > oo, 

trace,, h-E n trace,, h -> N |o, ^ J \^ hS fj ■ 

The last statement means convergence in distribution of random variables to a normal law with 
indicated expectation and variance. As noted in |3], Section 7, the result can be restated in terms 
of convergence of random fields to a Gaussian field on S with free boundary conditions. 

1.6. Derivation of Theorem ll.2l We now show, using the variational approach due to Johansson 
[13], that the Gaussian convergence stated in Theorem 11.21 follows from a generalized version of 
Theorem ll.il which we now state. 

Fix a real-valued test function h and consider the perturbed potentials 

Qn-~Q--h. 

n 

We denote by u n the one-point function of the density of states P„ associated with the potential 
Q n . We write a n for the measure with density fi„ and v„ for the measure n(a n - a), i.e., 

(1.3) 9 n [f] = na n (f) - no(f) = E„ Tr„[/] - no(f). 

Theorem 1.3. For all f e q°(C) we have 

v n {f)-v„{f)^^§Vf-Vh s . 
A proof of Theorem |1.3| is given in Section |4] 
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Claim. Theorem 1 1.2 1 is a consequence of Theorem 1 1.31 

Proof. Denote X n = Tr„ h - E„ Tr„ h and write a„ = E„X„. By Theorem ll.31 

On — > a where a = — \ V/z s • V/z s . 

More generally, let A > be a parameter, and let E„ ,\ denote expectation corresponding to the 
potential Q - (Ah)/n. Write 

F„(A) := log E„e AX ", < A < 1. 
Since F' n (A) - E„ ; aX„, Theorem [T3] implies 

(1.4) F,',(A) = E„, A X„ -» Aa, 
and 

(1.5) logE„e x » =F„(1) = f F;(A)dA -> ^, as n ^ oo. 

Jo z 

Here we use the convexity of the functions F H , 

f;;(a) = e„, a x2 - (e„, a x„) 2 > o, 

which implies that the convergence in \IA\ is dominated: 

= f;,(0) < F' n (A) < f n (l). 

Replacing /j by f/z where t e 1, we get E n (e tx ") —> e f "I 2 as m — > oo, i.e., we have convergence of 
all moments of X„ to the moments of the normal N(0, a) distribution. It is well known that this 
implies convergence in distribution, viz. Theorem 11.21 follows . □ 

I. 7. Comments. 

(a) Related Work. The one-dimensional analog of the weighted RNM theory is the more well- 
known random Hermitian matrix theory, which was studied by Johansson in the important 
paper [13] . Indeed, Johansson obtains results not only random Hermitian matrix ensembles, but 
for more general (one-dimensional) j3-ensembles. The paper |13| was one of our main sources of 
inspiration for the present work. 

In |3l, it was shown that the convergence in theorems 11.11 and 11.21 holds for test functions 
supported in the interior of the droplet. See also |6|. In J3j, we also announced theorems 11.11 and 

II. 2l and proved several consequences of them, e.g. the convergence of Berezin measures, rooted 
at a point in the exterior of S, to a harmonic measure. 

Rider and Virag [15J proved theorems 11.11 and 1 1 .21 in the special case Q(z) = |z| 2 (the Ginibre 
ensemble). The paper |8| contains results in this direction for |3-Ginibre ensembles for some special 
values of jS. 

Our main technique, the method of Ward identities, is common practice in field theories. In 
this method, one uses reparametrization invariance of the partition function to deduce exact 
relations satisfied by the joint intensity functions of the ensemble. In particular, the method was 
applied on the physical level by Wiegmann, Zabrodin et al. to study RNM ensembles as well 
as more general OCP ensembles. See e.g. the papers ll9ll , [20J, [21J, [22]. A one-dimensional 
version of Ward's identity was also used by Johansson in 1 13]. 
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Finally, we wish to mention that one of the topics in this paper, the behaviour of fluctuations 
near the boundary, is analyzed from another perspective in the forthcoming paper [4j. 

(b) Assumptions on the potential. We here comment on the assumptions (A1)-(A3) which we require 
of the potential Q. 

The C Cl ' assumption (Al) is natural for the study of fluctuation properties near the boundary 
of the droplet. (For test functions supported in the interior, one can do with less regularity.) 

Using Sakai's theory |18|, it can be shown that conditions (Al) and (A2) imply that dS is a 
union of finitely many C a ' curves with a finite number of singularities of known types. It is not 
difficult to complete a proof using arguments from [11J, Section 4. 

We rule out singularities by the regularity assumption in (A3). What happens in the presence 
of singularities is probably an interesting topic, which we have not approached. 

Without singularities the boundary of the droplet is a union of finitely many C"' Jordan curves. 
Assumption (A3) means that we only consider the case of a single boundary component. Our 
methods extend without difficulty to the case of a multiply connected droplet. The disconnected 
case requires further analysis, and is not considered in this paper. 

(c) Droplets and potential theory. We here state the properties of droplets that will be needed for 
our analysis. Proofs for these properties can be found in |16| and IflCfl . 

We will write Q for the maximal subharmonic function < Q which grows as log |z| + O(l) when 
|z| — > oo. We have that Q = Q on S while Q is C 1,1 -smooth on C and 

fi(z) = Q s (z) + G(2,oo) / zeC\S, 

where G is the classical Green's function of C \ S. In particular, if 

U°(z) = J log— L^Q 

denotes the logarithmic potential of the equilibrium measure, then 

(1.6) Q+U° = const. 

The following proposition sums up some basic properties of the droplet and the function Q. 

Proposition 1.4. Suppose Q satisfies (A1)-(A3). Then dS is a C a ' Jordan curve, Q e W 2/0O (C), and 
therefore 

dQ = (dQ) s . 

Furthermore, we have 

(1.7) Q(z) - Q(z) x 5(z) 2 , ztS, <5(z)^0, 
where b{z) denotes the distance from z to the droplet. 

(d) Joint intensities. We will occasionally use the intensity fc-point function of the process W n . This 
is the function defined by 

(k) P„(n- = i{^„nD(z y , £ )*0j) k 

*> = SS " - det(K„(z !V z ; )). , =1 . 

In particular, = nu n . 
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(e) Organization of the paper. We will derive the following statement which combines theorems ll.il 
and 11.31 (whence , by Lemma [l.6i it implies Theorem II. 2 ) . 

Main formula: Let v„ be the measure defined in dl-3| . Then 



(1.8) 



lim v„(f) = — 
«-»oo tin 



J^Af + J^fAL + J fN(L s ) ds + V/ s ■ Vh s 



Our proof of this formula is based on the limit form of Ward's identities which we discuss in the 
next section. To justify this limit form we need to estimate certain error terms; this is done in 
Section [3] In the proof, we refer to some basic estimates of polynomial Bergman kernels, which 
we collect the appendix. The proof of the main theorem is completed in SectionH] 

2. Ward identities 

2.1. Exact identities. For a suitable function v on C we define a random variable W+[c] on the 
probability space (C", P„) by 



j*k 



1 yi v(Aj) - v(A k ) 

2 lr* A; - h 



In Tr n [vdQ] + Tr„[dv]. 



Proposition 2.1. Let v : C — > C be Lipschitz continuous with compact support. Then 

E„W„>] = 0. 

Proof. We write 



W^[v] = I n [v]-II n [v] + III n [v] 



where (almost everywhere) 

" t>(zy) -v(z k ) 



i«b>Uz) = \Y J 



Zj - Zfc 



; II n [v](z) = 2Y j dQ(z j )v(z j ) ; III n [v](z) = £ dv (zy) . 



j*k I * /=1 ;=1 

Let e be a real parameter and put Zy = (p (ij) = Cy + £-V (Cy) /2, 1 < ;' < n. Then 

d z Zj = ( | ty (Cy) \ 2 + \d(p (Cy) | 2 ) d% = [ 1 + e Redv (Cy) + O (e 2 ) ] d%, 



so that (with III„ = III„[v]) 



dA n (z) = [ 1 + e Re ffl„ (Q + o(e 2 ) ] dA„(Q. 



Moreover, 



log z; - zy = log C; - Cy + log 



= log I Ci - Cy I + £ Re 



e i>(C0-?(Cy) 



Ci-Cj 



0(e% 



so that 



(2.1) ^ log | z ; - - Zfc | _1 = ^log|C / -0c|~ 1 -£Re J„(Q + 0(e 2 ), as e -> 0. 

7** /** 
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Finally, 

Q(z,) = Q(C; + ^(C;)) = Q(C;) + £ Re (<?q(q) *(C/)), 

so 

n n 

(2.2) 2« Yj Q ( z /) = 2 " E Q fe) + £ Re + ° H ■ 

7=1 ;=i 

Now l|2.H l and l|2,2|l imply that the HamiltonianH„(z) = log | z ; - - z^ | 1 +2« £J =1 Q(z ; ) satisfies 

(2.3) H n (z) = H n (Q + e • Re (-I n (Q + II n (Q) + O (e 2 ) . 
It follows that 

Z„ : = J e" H »< z > dA„(z) = £ e" H »< c >- £ Re H-(0 + i4(Q) + o(^) [ a + £Reffl(i( Q + o(e 2 )] dA„(Q. 

Since the integral is independent of e, the coefficient of e in the right hand side must vanish, 
which means that 

(2.4) Re f (III n (Q + I n (Q - II n (Q) e" H » (C) dA n (Q = 0, 

JC" 

or Re E„ W,~J"[c] = 0. Replacing v by ic in the preceding argument gives Im E„ W,~J"[c] = and the 
proposition follows. □ 

Applying Proposition 12. II to the potential Q„ = Q - h/n, we get the identity 

(2.5) E„W„>] = 0, 
where 

(2.6) W+[v] = W„ + M + 2 Ucwbdfe]. 
If we denote 

we can rewrite (12.51 1 and 02. 6t as follows, 

(2.7) E n B„M = 2E„ Tr n [vdQ] - a n (dv + 2vdh), 
where we recall that &„ is the measure with density u n . 

2.2. Cauchy kernels. For each z e C let k z denote the function 

k z (A) ' 



z-A' 

sozh a{k z ) is the Cauchy transform of the the measure a. We have (see (|1.6t ) 

tx(fc z ) = 2«9Q(z). 

We will also consider the Cauchy integrals a n (k z ) and a n (k z ). We have 

<?z[0n(fc 2 )] = TlU n (z), d z [0 n (k z )] = TTM„(z), Z € C, 

and 

a„(k z ) -> a(k z ) 
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with uniform convergence on C (the uniform convergence follows easily from the one-point 
function estimates in Lemma |3H and Theorem l3.2t . 

Let us now introduce the functions 

D„(z) = v n {k z ) ; D„(z) = v n (k z ). 

We have 

(2.8) D„(z) = n[a n (k z ) - 2dQ(z)] , dD n = nn(u„ - u), 
and if f is a test function, then 

(2.9) v n {f) = \ j fdD n = -i f 5f-D n . 

Let K„ denote the correlation kernel with respect to Q„. Using D„, we can rewrite the B n [v] term 
in the Ward identity as follows. 

Lemma 2.2. One has that 

E„B n [v] = 2 J vdQ-K n + J vD n u n -^ IT V -^-^^ \K n {z,w)\ 2 . 

(In the first integral K„(z) means the 1-point intensity R„ (z) = K„(z,z).) 



Proof. We have 



where 



2n JJ C2 



R^\z,w) = K„(z)K„(ro) - \K n (z,zv)\ 2 . 



The integral involving K n (z)K n (w) is 

i rr _^_ Kn(z)Rn(zt7) = f ^ . k„( Z ) • a„(/c z ), 

" J JC2 z - w Jc 
and by gS> a n (k z ) = \t) n + 25Q. 

2.3. Limit form of Ward's identity. The main formula | |1.8t will be derived from Theorem 
below. In this theorem we make the following assumptions on the vector field v: 

(i) v is compactly supported; 

(ii) v is Lip-continuous in C; 

(iii) v is uniformly C 2 -smooth in C \ dS. 

(The last condition means that the restriction of V to S and the restriction to (C \ S) U dS are both 
C 2 -smooth.) 

Theorem 2.3. Ifv satisfies (i)-(iii), then asn — > oo, 

- f vD„d5Q + - f v(dQ-dQ)3D n ^-\a{dv)-2o{vdh). 
n Js 71 Jc\s 2 
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Before we come to the proof, we check that it is possible to integrate by parts in the second 
integral in Theorem l2.3l To control the boundary term we can use the next lemma. 



Lemma 2.4. For every fixed n we have 



Proof. We have 



Since 



D„(z) 



/(u n - u)d z A _ C 
z-A = J 



1 



1 1 



z - A 
A 



(u„ - u)d A 



z - A z z 2 1 - A/z' 



we need to show that the integrals 



|A| |w„ - u| 



d 2 A 



U-A/z| 

are uniformly bounded. To prove this, we only need the estimate u„ (A) < jjp which holds (for suf- 
ficiently large n) by the growth assumption l|1.2t and the simple estimate u„(A) < C exp(-2n(Q(A) - 
Q(A))), which is given below in Lemma [3. II □ 

Using that dQ - dQ in the interior of S, we deduce the following corollary of Theorem |2.3l 

Corollary 2.5. ("Limit Ward identity") Suppose that v satisfies conditions (i)-(iii). Then as n — > oo we 
have the convergence 

| [vddQ + 5v(dQ - dQ)] D„ ^ ~a(dv) - 2o(vdh). 

2.4. Error terms and the proof of Theorem 12.31 Theorem 12 . 3 1 folio ws if we combine the expres- 
sions for £„£>„[£>] in l|2.7t and Lemma 12.21 and use the following approximations of the last two 
terms in Lemma [2. 21 More precisely, if we introduce the first error term by 

(2.10) - [f p < 2 )- p < ro ) \K n (z,wf = a n (dv) + ellv], 

n J J z-w 

and the second error term by 

(2.11) £>] = nj vD n (u n - u) =~\fdv ^. 
Using \2.7\ , Lemma |Z2] and that dQ = dQ a.e. on S, one deduces that 

(2.12) - fvD n ddQ+- f v{dQ-dQ)dD n = -\o{dv)-2o{vdh)+\e l n [v}--el[v}+o{\), 

where o(l) = (a- d n )(dv/2 + 2vdh) converges to zero as n — > oo by the one-point function estimates 
in Lemma [3. II and Theorem l3.2l 

In the next section we will show that for each v satisfying conditions (i)-(iii), the error terms 
e'nlv] tend to zero as n — > oo, which will finish the proof of Theorem |2,3l 
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3. Estimates of the error terms 

3.1. Estimates of the kernel K n . We will use two different estimates for the correlation kernel, 
one to handle the interior and another for the exterior of the droplet. 

(a) Exterior estimate. Recall that K„(z, w) is the kernel of the n-point process associated with 
potential Q„ = Q - h/n; as usual, we write K„(z) = K„(z,z). We have the following global 
estimate, which is particularly useful in the exterior of the droplet. 

Lemma 3.1. For allz e C we have 

K„(z) < ne- 2n ®-®<*\ 
where the constant is independent ofn and z. 



This estimate has been recorded (see e.g. [[2], Section 3) for the kernels K„, i.e. in the case h = 0. 
Since obviously 

Wp\\ e -2nQ„ - Wp\\e-2nQ, 

we have K„(z) x K„(z) with a constant independent of z. Indeed, K„(z) is the supremum of 

i i2 

p(z) e~ 2 "Q( z ) where p is an analytic polynomial of degree less than n such that ||p|| e -2nQ < 1, and we 
have an analogous supremum characterization of K„(z). Hence the case h + does not require 
any special treatment. 

In the following we write 

S(z) = dist(z,<?S) 

and 

log 2 n 

&n = — 

yn 

By our assumption on the droplet (see Proposition [TT4J we have 

Q(z) - Q(z) > 6 2 (z), z $ S, 6(z) -> 0. 

In view of the growth assumption I jl.H , it follows that for any N > there exists Ov such that 
K„(z) < Cn«~ n when z is outside the 5 n -neighborhood of S. 

(W Interior estimate. Recall that we assume that Q is real analytic in some neighbourhood of 
S. This means that we can extend Q to a complex analytic function of two variables in some 
neighbourhood in C 2 of the anti-diagonal 

{(z,z) : zeS}cC 2 . 

We will use the same letter Q for this extension, so 

Q(z) = Q(z,z). 

We have 



Q(z,w) = Q(w,z) 

and 

dxQ(z,z) = dQ(z), d l d 2 Q{z,z) = ddQ{z), d\Q{z,z) = d 2 Q(z), etc. 
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With the help of this extension, one can show that the leading contribution to the kernel K„ is of 
the form 

(3.1) K(z,w) =-{d l d 2 Q){z,w)n^ 2Q{z ' m) - Q(z) - Q(w) \ 

TL 

In particular, we have 

2tl 

We shall use the following estimate in the interior. 
Theorem 3.2. Ifz e S, 5(z) > 2b n , and if\z - w\ < <5„, then 

\K n (z,w)\ = |k#(z,m;)| + 0(1), 
where the constant in O(l) depend on Q and h but not on n. 

Similar types of expansions are discussed e.g. in |5], |1J, [2]. As there is no convenient reference 
for this particular result, and to make the paper self contained, we include a proof in the appendix. 

We now turn to the proof that the error terms £„[t/j and £ 2 [v] are negligible. See | |2.10t and I j2.11^ . 
Our proof uses only the estimates of the kernels K„ mentioned above. Since the form of these 
estimates is the same for all perturbation functions h, we can without loss of generality set h = 0, 
which will simplify our notation - no need to put tildes on numerous letters. 

3.2. First error term. We start with the observation that if w e S and 5(w) > 25„ then at short 
distances the so called Berezin kernel rooted at w 

(w) _ \K n (z,w)\ 2 
" 1 ' K n (zv,w) 

is close to the heat kernel 

H ( ™\z) = -cne- aiz - w ? , c := 2ddQ(w). 

TL 

Both kernels determine probability measures indexed by w. Most of the heat kernel measure is 
concentrated in the disc D(w, 5„), 

f H<;"\z)dA(z) < 4, 

JC\D(H),6„) n 

where N denotes an arbitrary (large) positive number. 

Lemma 3.3. Suppose that w e S, 5(w) > 25 n and \z — w\ < b n . Then 

\Bt"\z)-H^\z)\ < nb n . 

Proof. By Theorem |3.2| we have 

B M( Z) = %^ + (1). 
K„(if, w) 
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Next, we fix w and apply Taylor's formula to the function z i— > K # (z, w) at z = w. Using the explicit 
formula (I3.lt for this function, and that 



Q(z,w) + Q(z,w) - Q(z,z) - Q(w,w) = [Q(z,w) - Q(w,w)\ + [Q(w,z) - Q(z,z)] = 
9Q(w)(z -w) + ^d 2 Q(w)(z - wf + dQ(z)(w - z) + \d 2 Q{z){z -wf + --- = 
[dQ{w) - dQ(z)](z -w) + d 2 Q{w){z -wf + --- = -ddQ{w)\z - w\ z + . . . 



we get 



|K*(z,i(;)| 2 1 



K^(w,w) TL 



= H(C + 0(|Z - W\)] ne -"'\z-™\ 2 +0(n\z-z»\ 3 ) 

= H^\z) + 0(n\z - w\), 



and the assertion follows. 

Corollary 3.4. If w e S and 5(w) > 25 n , then 

B< w >(z) dA(z) < nb\ = o(l). 



f 



JC\D(wA,) 

Proof. We write D„ = D(w, 5„) and notice that 



r si a,> = i - r b< w) = i - f + r (H< n,> - b<t>) = r + r (H< a,> - b<t>). 

Jc\d„ Jd„ Jd„ Jd„ Jc\d„ Jd„ 



The statement now follows from Lemma [3. 3 1 

Proposition 3.5. Ifv is uniformly Lipschitz continuous on C, then e\[v] — » as n — » oo. 
Proof. We represent the error term as follows: 



J fa] = I «„(w)F n (z0) ; F„(io) = I 

J(a>) J(z) L 



v(z) - v(w) 
z — w 



dv(w) 



Bf\z). 



By the assumption that v is globally Lipschitzian, we have that 



f 

J 6 



u„(h;)F r (io) 



(6(h?)<2S„) 



< I u„(w) = 0(1). 

'{S(a>)<26„) 



If 6(h;) > 2S„, then 



I 



zeD(z»,5„) 



d(z) - v(w) 



z — w 



- <9c(a;) 



B<r>(z) 



+ const. 



/ 



Bt\z), 



where the last term is o(l) by Corollary |3.4| Meanwhile, the integral over D(w, 5„) is bounded by 



IX 



z£D(w,6 n ) 



v(z) - v{w) 
z-w 



- dv{w) 



H^(z) 



J DC 



-const. | \Bt\z)-H^\z)\, 

D(w,6„) 



where we can neglect the second term (see Lemma l3.3|l . Finally, 

viz) — v(w) - z — w 

- dviw) = dv(w) + o(l), 

z-w z-w 

(this is where we use the assumption v e C J (S)), so the bound of the first term is o(l) by the radial 
symmetry of the heat kernel. □ 
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3.3. Second error term. We shall prove the following proposition. 
Proposition 3.6. Ifv is uniformly Lipschitzian, then 

i r - d 2 

The proof will involve certain estimates of the function 

Jc Z -C 

(Here K # (Q = nu ■ Is.) It will be convenient to split the integral into two parts: 
D„(z) = C(z) + K„(z) := (JT + ) K " (C ^^" (Q d 2 C, 

where 

B„ = {z : 6(z) < 26„}. 
By Theorem l3.2l and Lemma [3.1l we have |K„ - K*| < 1 in C \ B„, and therefore 

|R„I < 1. 

Hence we only need to estimate C n , the Cauchy transform of a real measure supported in B„ - a 
narrow "ring" around <9S. We start with a simple uniform bound. 

Lemma 3.7. The following estimate holds, 

\\Dnh- £ V«log 3 «. 

Proof. This follows from the trivial bound |K„ - K*| < n and the following estimate of the integral 

r JK_ 

Jb„ iz-cr 

Without loosing generality we can assume that z = and replace B n by the rectangle \x\ < 1, 
|y| < 6„. We have 



Jb„ ici j-i j-6„ v* 2 + v 2 



V* 2 + y 2 

where I is the integral over {|x| < 5„] and is the remaining term. Passing to polar coordinates 
we get 

,<5„ 



r n rdr 
Jo r 



and 

r 1 dx 

II<5 n — x<5„|logS„|. 



Lemma |3~7| gives us the following estimate of the second error term, 
(3.2) |£>]|<log 6 n|Mhp, 
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which comes rather close but is still weaker than what we want. Our strategy will be to use ^ .._ , 
and iterate the argument with Ward's identity. This will give a better estimate in the interior of 
the droplet. 

Lemma 3.8. We have that 



\C„(z)\ < 



log 6 M 

<5(z) 3 ' 



z e S. 



we 



Proof. Let ip be a function of Lipschitz norm less than 1 supported inside the droplet. Then 
have 

eltylSl, tfji/^logS, 
where the constants don't depend on ip. (The first estimate follows from Proposition |3.5i and the 
second one is just 1 13.21 ). This means that the error e n := + e 2 in the identity l|2.12|l is bounded 
by log 6 n for all such ip, i.e. (since dQ = dQ a.e. on S), 



xpD n AQ 



: 0(1) + |e„[^]| <log 6 n, 



and therefore, 



ipC n AQ 



: log 6 n 



For zeS„, we now set 25 = b(z) and consider the function 

fs-ie-z 



ip(Q = max 



AQ(Q 



, 



Then ip has Lipschitz norm « 1, and by analyticity of C„ we have the mean value identity 

J \pC n AQ = 2ttC„(z) J^(6 - r)rdr = nC n (z)5 3 /3. 
We conclude that \C„(z)\ < 5~ 3 log 6 n. 



Finally we need an estimate of C„ in the exterior of the droplet. This will be done in the next 
subsection by reflecting the previous interior estimate in the curve T :- dS. 

Let us fix some sufficiently small positive number, e.g. £ = j$ will do, and define 

y„ = n~ £ . 

Denote 

r n = {C + 7„v(Q: CeT}, 
where v(Q is the unit normal vector to T at C, e T pointing outside from S. We will write intr„ 
and extr„ for the respective components of C \ T n . In the following, the notation a < b will mean 
inequality up to a multiplicative constant factor times some power of log n (thus e.g. 1 < log 2 n). 

Let L 2 (r„) be the usual I 2 space of functions on T n with respect to arclength measure. We will 
use the following lemma. 

Lemma 3.9. Suppose that v is uniformly Lipschitzian. Then 

Given this estimate, we can complete the proof of Proposition |3.6| as follows. 
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Proof of Proposition ^. 61 Applying Green's formula to the expression for e%[v] (see l|2.11t ), we find 
that 

S -\\D n \\ 2 LW ) + -H D "HL 2 (r„)' 
where intr„ is the y n -neighbourhood of the droplet. The second term is taken care of by Lemma 
13.91 To estimate the first term denote 

A„ = {6(z) < y„). 

The area of A n is x y ni and in S \ A n we have |D„(z)| < y~ 3 (Lemma 13. 8t . We now apply the 
uniform bound |D„| < yfn in A n (Lemma 13. 7) . It follows that 

H D "llL2(intr„) = f + \ <n\A n \ + y- n 6 , 
Ja„ Js\a„ 

whence 

ll D »llL 2 (intT„) = °(»)- 

This finishes the proof of the proposition. □ 
3.4. Proof of Lemma [3^91 Let us first establish the following fact: 
(3.3) |lm[v(C)C„(C + 7nV(Q)]| < <Wn, CeT. 



Proof. Without loss of generality assume that C = and v(Q = i. 

The tangent to T at is horizontal, so T is the graph of y = y{x) where y{x) = 0(x 2 ) as x — > 0. 
We will show that 

(3.4) |Re[C„(z>„) - C„(-iy„)]\ < yfny^. 

This implies the desired estimate d3 .3t , because by Lemma |3~8l 

\C„{-iy„)\ < y~ 3 < yfny^. 

To prove (13. 4t we notice that 

1 1 



I:=Re[C n (iy n )-C n (-iy n )]= f Re 

Jb„ 



z - iy„ z + iy n 



p„(z)dA(z), 



where we have put p„ = K„ - K*, viz. \p n \ < n. 

We next subdivide the belt B n = {6(z) < 25,,} into two parts: 

B' n = B n n {|x| < Vr«), K = B n \ B' n . 

Clearly 



(3.5) 



\I\< 



Jbl 



1 



1 



z - iy n z - iy n 



z - zy„ z + zy„ 



The integral over B' n in the right hand side of < |3 - 5b is estimated by 

Jb;, Jb;, * 2 + yf, 
because if z = x + z'y 6 B', then |y| < x 2 + 6„ < x 2 + y„. 
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We estimate the integral over B" in d3-5b by 

r x f 5 < 6 «> / " r $ x 6 « vy«- 

Jb;; Jb;; x ^VF» 

It follows that 

|J| -< n<5„ Vy« < y/ny7i- 

This establishes [ 13.4t , and, as a consequence, I l3,3t , □ 

To finish the proof of Lemma [3. 9 1 we denote by v„(-) the outer unit normal of Y n . Using < |3.3b and 
Lemma |3~7l we deduce that 

|Im v„C„| < yny„ on Y„. 
Next let D» be the exterior of the closed unit disk and consider the conformal map 

(p„ : ext(r„) — > D„ oo i— > oo. 

We put 

_ (p„C„ 
- T, ■ 



Then F„ is analytic in ext(r„) including infinity, and we have 

|2 



(3.6) ||ImF„||j 2 < y[ny^. 



To see this note that 

Im[v„C„] 



Im F„ 



and recall that we have assumed that Y is regular (A3), which means that \(p' n \ is bounded below 
by a positive constant. 

Now note that (p n (z) / (p' n (z) = r n z + O(l) as z — > oo, where the r n are uniformly bounded. This 
gives 

F„M = r n f (Kn-K*) = r„ f (K* - K„) = O(l), 

Jb„ Jc\b„ 

where we have used Lemma |3. II and Theorem 13.21 to bound the integrand. Therefore, by 
since the harmonic conjugation operator is bounded on the Hardy space H 2 , 

||ReF„||^ 2(rn) x ||ImF„||^ 2(rn) + O(l) < yfny^. 

This completes the proof of Lemma [379] 



4. Proof of the main formula 

In this section we will use the limit form of Ward's identity (Corollary 12. 5t to derive our main 
formula (|1.8t : for every test function / the limit v(f) := lim^oo v n (f) exists and equals 



(4.D Hf) = I 



J A/ + JfAL + fN(L s )ds + V/ s ■ V/; s . 
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4.1. Decomposition of f. The following statement uses our assumption that dS is a (real analytic) 
Jordan curve. 

Lemma 4.1. Let f e C^°(C). Then f has the following representation: 

f = f ++ f_ + f 0/ 

where 

(i) all three functions are smooth on C; 

(ii) 9f + = and df. =0inC\ S; 
(Hi) f ± = O(l) at oo; 

(iv) fo = on dS. 

Proof. Consider the inverse conformal maps 

(p : D» -> C \ S, ip : C \ S D», co ^ co, 
where E), = {|z| > 1}. On the unit circle T, we have 

oo 

F:= /o^^„C» € r(T). 

— oo 

The functions 

oo 

F + (z) = £ a n z", F_(z) = £ J, (z e D.) 

— oo 1 

are C°° up to the boundary so we can extend them to some smooth functions F+ in C. The 
conformal map ip also extends to a smooth function ip : C — * C. It follows that 

f ± := F ± o xp e q°(C), 

and /+ satisfy (ii)-(iii). Finally we set 

/o=/-/ + -/_. 

□ 

Conclusion. It is enough to prove the main formula (|4.1t only for functions of the form / = 
/++/_+ fo as in the last lemma with an additional assumption that fo is supported inside any 
given neighborhood of the droplet S. 

Indeed, either side of the formula l|4.1t will not change if we "kill" fo outside the neighborhood. 
The justification is immediate by Lemma [3711 

In what follows we will choose a neighborhood O of S such that the potential Q is real analytic, 
strictly subharmonic in O, and 

dQ + dQ in 0\S, 

and will assume supp(/ ) c O. 
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4.2. The choice of the vector field in Ward's identity. We will now compute the limit 

v(/):=limv„(/) 
(and prove its existence) in the case where 

/ = / + +/o- 

To apply the limit Ward identity 

(4.2) ^fjvddQ + dv(dQ-d&]D n ^~o(dv)-2o(vdh), (n ^ oo), 

(see Corollary 12. 5 1 , we set 

V = V+ + V , 



where 



and 



This gives 



But in C \ dS we have 



dfo 1 fo 1 

<?<?Q <?Q - <?Q 

C+ = • lc. 

ddQ 



df fo 
v = ■ l s + — J - }! —r ■ l C \s- 



vddQ + do-d(Q-Q) = df, 
so comparing with (12. 9t , we find that 

(4.3) - 2v n (D = \ f H<5Q + dv(dQ ~ dQ)] D„. 

However, to justify that l|4.2t holds, we must check that v satisfies the conditions (i)-(iii) of 
Corollary |2.5l 

Lemma 4.2. The vector field v defined above is Lip(<C) and the restrictions ofvtoS and toS, : = (C \ S) U dS 
areC™. 



Proof. We need to check the following items: 

(i) v\s, is smooth, and (i') v\$ is smooth; 

(ii) Vq is continuous on dS, and (ii') same for v+. 

The items (i') and (ii') are of course trivial. (E.g., 5f + ■ 1$ = df + ) 

Proof of (i). We have v = fo/g in C \ S where g = dQ - dQ. Since the statement is local, we consider 
a conformal map <p that takes a neighbourhood of a boundary point in S onto a neighbourhood 
of a point in ]R and takes (parts of) dS to ]R. If we denote F = fo ° <p and G = g o <fi, then F = and 
G = on ]R. Moreover, G is real analytic with non-vanishing derivative G y . Thus it is enough to 
check that 

H(x, y) = F -^-lA 

y 
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has bounded derivatives of all orders. We will go through the details for H, H y , H yy , Applying 

the same argument to H x , we get the boundedness of the derivatives H x , H xy , H yy , . . . , etc. 

Let us show, e.g., that H' :- H y is bounded. We have 

tji yF'-F y(F' + O(y))-(yF' Q + O(if)) 

H = — = j = °®. 

where F' :- F'(-, 0) and all big O's are uniform in x. (They come from the bounds for the 
derivatives of F.) Similarly, 

y 2 F" - 2yF' + 2F 



H" = 

y 3 



The numerator is 



y 2 (F^' + 0(y)) - 2y(F' Q + yF'J + 0(y 2 )) + 2(yF' Q + iy 2 F[,' + 0(y 3 )) = 0(y 3 ), 
etc. (We can actually stop here because we only need C 2 smoothness to apply Theorem |2.3l ) 

Proof of (ii). Let n = n(Q be the exterior unit normal with respect to S. We have 

/o(C + 6n)~65 n /o(0 = 26-(5/o)(C)-KO, as 6 | 0. 
Similarly, if g := dQ - dQ, so g = on dS and dg = ddQ in C \ S, then 

g{t, + bn)~bd n g{Q=2b-{dg){Q-l^Q, as 5 1 0. 

It follows that 

/o(C + 6») 5f (Q 
g{t + bn) ~d5Q(Q' [ 1 h 
which proves the continuity of Vq. □ 

We have established that v = vo + v+ satisfies conditions (i)-(iii) of Corollary 12.51 Thus the 
convergence in l|4.2t holds, and by [ 14.31 1 we conclude the following result. 



Corollary 4.3. Iff = f + /+, then 



1 

v(f) = -a(dv) + a{vdh). 



4.3. Conclusion of the proof. 

(a). Let us now consider the general case 

f = f ++fo+ f_. 

By the last corollary we have 

l 5f+ 
Hf+) = ^o(dv+) + a(v + dh), v + := ^- ■ l s . 

Using complex conjugation we get a similar expression for v(/_): 

1 - - df- 

v(f_ ) = -oidvJ) + o(v-dh), v- ■- ■ l s . 
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1 1 . 

v(/_) = {*(/_) = -a(dv-) + a(v-dh) = -a(dv-) + a{v-dh). 
(Recall that h is real-valued.) 
Summing up we get 

(4.4) v(f) = \ [a{dv ) + a{dv + ) + o{dv-)\ 
and 

(4.5) v(f) - v(f) = a{v + dh) + a(v dh) + o(v-dh). 



(b) Computation ofv(f). Recall that 



da{z) = — AQ(z)l s (z)cfz , L = log AQ. 

2 71 



Using l|4.4t we compute 



H3Q 



--ll d5f -ll 5f ° dL -hl 5f + dL -il df - 5L 



ddQ -~) 2ti J s A?Q " ^ ' ~' 2ti J s " \^Q/ " - ~ 2ti Js [ddQj 



ddQ 



At this point, let us modify L outside some neighborhood of S to get a smooth function with 
compact support. We will still use the notation L for the modified function. The last expression 
clearly does not change as a result of this modification. We can now transform the integrals 
involving L as follows: 

- f df dh- f df + dh- f df- 3L= ffo ddL + f(f++ /-) <?<9L 
Js Jc Jc Js Jc 

= \fddL+ \ f ddL, 
Js Jc\s 



and we conclude that 



v(/) = 



8tt 



f A/ + f /AL + f f s AL 
Js Js Jc\s 



Note. The formula for v(f) was stated in this form in 1 3 1 . 
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f f s AL= f (f s AL-LAf s ) 

Jc\s Jc\s 

= [ if- d n L - d n f ■ L S ) ds 

JdS 

= f (f s -d n L-f s -d„L s ) ds 

JdS 

= f fN(L s )ds 

JdS 



In conclusion, 
(4.6) 



•Af) 



8tt 



f A/ + f/AL + f fN(L s 
Js JS JdS 



(c) Computation of[v(f) - v(f)]. Using the identity (14.5t , we can deduce that 



v(/)-v(/) = - 



71 



f 3f + dh+ f dfSh+ f 5f dh 
\J s *j s *j s 



JS JS 

fS . X7hS 



This is because 



and similarly 



On the other hand, 



Therefore, 



and this is equal to 



\ V/ ■ V/z + | V/ s ■ V/z , 

Js Jc\s 



v(/)-v(/) = 



— r v/ s -v/2=— r vf s -v/z s . 



Applying (14. 6t we find that 

1 



'(/) 



8tt 



A/ + J^/AL + £ /7V(L S ) + V/ s • Wz s , 



and the main formula fl4.lt has been completely established, q.e.d. 
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Appendix: Bulk asymptotics for the correlation kernel 

Polynomial Bergman spaces. For a suitable (extended) real valued function cp, we denote by 
L 2 the space normed by ||/|| 2 = J £ \f\ 2 e~ 2c ^. We denote by A 2 the subspace of L 2 consisting of a.e. 

entire functions; !P„(e~ 2<f) ) denotes the subspace consisting of analytic polynomials of degree at 
most n-1. 

Now consider a potential Q, real analytic and strictly subharmonic in some neighborhood of 
the droplet S, and subject to the usual growth condition. We put 

Q = Qn = Q- -h, 
n 

where h is a smooth bounded real function. 

We denote by K the reproducing kernel for the space P n {e~ 2nQ ), and write K w (z) = K(z, w). The 
corresponding orthogonal projection is denoted by 

P n :L 2 nQ ^p n (e- 2 "Q) : f»(f,K w ) nQ . 

The map P n : L 2 ^ — > r P n {e~ 2n ®) is defined similarly using the reproducing kernel K for the space 
Pn(e- 2nQ )- 

We define approximate kernels and Bergman projection as follows. In the case h = 0, the well- 
known first order approximation inside the droplet is given by the expression 

Kl{z) = -{d 1 d 2 Q){z,w) ne 2n Q^, 
n 

where Q(-, •) is the complex analytic function of two variables satisfying 

Q{w,w) = Q{w). 

If the perturbation h + is a real-analytic function, we can just replace Q by Q in this expression. 
Note that in this case, the analytic extension h(-, •) satisfies 

h(z, iv) = h(w) + (z - w)dh(w) + (z — > w). 

This motivates the definition of the approximate Bergman kernel in the case where h is only a 
smooth function: we set 

g*(z) = K* w (z) e- 2h -^\ 

where 

h w (z) := h(w) + (z - w)dh(w). 
The approximate Bergman projection is defined accordingly: 

P*J(w) = (f,K* w ) nQ . 



The kernels £^(z, w) do not have the Hermitian property. The important fact is that they are 
analytic in z. 
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Proof of Theorem |3.2l We shall prove the following estimate. 
Lemma 4.4. Ifz e S, S(z) > 2b„, and if\z - w\ < 5 n , then 

\K w (z) - Kt(z)\ < e' lQ{z) e nQ(w) . 
Before we prove the lemma, we use it to conclude the proof of Theorem |3.2| Recall that 

K„(z,h>) = K w (z) e~ nQ(z) e~ nQiw) . 

If we define 

¥*{z,w) = X*(z) e""Q (z) e - nQ{m \ 

then by Lemma 23) 

K n (z,w) = K*(z,H;) + 0(l). 

On the other hand, we have 

K # n (z,iv) = Kl(z,w) ^)+hdo)-2h-M) ^ 

so 

|«*(z,tr)| = KM\ (1 + 0(\w - z\ 2 ) = \K # n (z,w)\ + 0(1). 

It follows that 

\K n (z,w)\ = \K # n (z,w)\ + 0(l), 

as claimed in Theorem |3.2| □ 
It remains to prove Lemma [4.41 

Lemma 4.5. Iff is analytic and bounded in D(z; 25„) then 

\f(z)-Pt( Xz f)(z)\<^-e'^\\f\\ nQ . 

y/n 

Here X — Xz is a cut-off function with % — 1 m C(z;36„/2) and j = outside D(z;26 H ) satisfying 

115*112*1. 

Proof. Wlog, z = 0, so Pf,(xf)(0) is the integral 

f = I J x( q . /(Q . 2(^0(0, ■ e 2WO-K0)-C5A(0)] . ne -2 n [Q(CC)-Q(o,C)]_ 

Since 

5c j e -2„[ Q( a)-Q(o,C)]] = _ 2 [,9 2 Q(C, Q - <? 2 Q(0, 0] TO -2»I«CO-Q(P / 0] / 
we can rewrite the expression as follows: 



= ~f\ /(CMCM(C)B(Q 5 [ e -2»IQ(C-0-Q(o,C)]] / 



where 

A(Q= ^ 5 i^Q)(0,Q 



d 2 Q(C,Q-d 2 Q(0,Q' 
and 

B(Q = e 2[ft(C)-M0)-C*(0)]_ 
A trivial but important observation is that 

A,B = 0{1) , 5A = 0(\Q) , 5B = 0(\Q), 
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where the O-constants have uniform bounds throughout. 
Integrating by parts we get 

f = f(0)+e 1+ e 2 , 

where 

£l = J L ^ AB e -2«[Q(C,C)-Q(0,C)] / g 2 _ J f X ^ AB ) e -2n[Q(C)-Q(0,C)]_ 

Using that 

| £1 |<1 J | /M ^| e -2n[Q(Q-ReQ(0,C)] / | £2 |< J ^ | /( e -2n[Q(C)-Re 0(0,0^ 

and noting that Taylor's formula gives 

e -n[Q(Q-2ReQ(0,Q] < e «Q(0)-c„|C| 2 ^ ^ „ AQ(0) > 0), 

we find, by the Cauchy-Schwarz inequality (since |CI > <5« when dx(Q + 0) 

p -cn&l _ -i 

|ei|e-" Q(0) < -^\\f\\n Q \\dx\y < -pWfWnQ 
On yn 

and 



The proof is finished. 



Suppose now that dist(z, C \ S) > 2<5„ and \w -z\< 5„. 
From Lemma [4.51 we conclude that 

(4.7) \K w (z) - Pn [xzK] (z)| < e" Q(z) e"»"\ 

This is because 



^ [xzK z ] (ip) = Gfc&^nQ = te£W„Q = P» (z), 



so 



|^(z)-n[xzK#](z) = K 2 (i(;)-P#[x 2 K z ](i<;) 
and because (cf . [2], Section 3) 

lliy = < V^" Q(z) - 

On the other hand, we will prove that 



(4.8) 



£»-p« [rf](w) 



«Q(z) _nQ(a;) 



which combined with 1)4.71 proves Lemma [4.41 The verification of the last inequality is the same 
as in ||5| or (TJ, depending on the observation that L*q = with equivalence of norms. We give 
a detailed argument, for completeness. 

For given smooth /, consider u, the 1?^ -minimal solution to the problem 

(4.9) 5u = 3f and w-/ef»„_i- 

Since ||m||„q < ||u|| n Q, the L^-minimal solution u to the problem J4.9I I satisfies ||m||„q < C||u||„q. We 
next observe that P n f is related to the L^-minimal solution u to the problem | |4.9t by u = f - P n f. 
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We write 

u(Q = Xz (QKi(Q-P n [xzK # w ] (Q, 
i.e., u is the L^-minimal solution to (|4.9t for f = Xz' Let us verify that 

(4.10) IMI„q<||^(av£*)|L q . 

To prove this, we put 

2</)(Q = 2Q(Q + «- 1 log(l + |C| 2 ), 

and consider the function Vq, the L^-minimal solution to the problem 5v = d(xz • K^)- Notice 
that <p is strictly subharmonic on C. f3y Hormander's estimate (e.g. [12], p. 250) 

lk<< j c \d( Xz -kl)\ _ 



Since Xz is supported in S, we hence have 



1 1 1 - II 

We next observe that by the growth assumption on Q near infinity, we have an estimate 
ncf) < nQ + const, on C, which gives ||i>oI|«q < ||i>olln</>- It yields that 



M«Q<^||5(xz<)l • 



But c - Xz ■ K-f„ belongs to the weighted Bergman space A 2 ,. Since 2n(p(Q = (n + 1) log |CI 2 + O(l) 
as C — * 00 / the latter space coincides with *P n -i as sets. This shows that vq solves the problem 
(14.91 . Since ||m||„q < IbolInQ/ we then obtain l|4.10t . 
By norms equivalence, J4.10t implies that 

(4.11) M\nQ^^r¥(xz-K)\ ^ 

yn II v '"nQ 

where 

U = X zKt-Pn [xzK] 

is the L^-minimal solution to | |4.9t with / = XzKf v . 
We now set out to prove the pointwise estimate 

(4.12) \u(z)\ < ne- cn6 "e n(Q(z)+Q(w) \ 
To prove this, we first observe that 

da(Q = 5(xz ■ K) (0 = dXz{Q ■ KiQr 

whence, by the form of K* and Taylor's formula, 

\3u{QU~ 2 " Q(Q £ n 2 \5xz(Q\ 2 e MQ{w) ~ clC ~ wl2) 



with a positive constant c ~ AQ(z). Since |C — w | > 5 n /2 when dx(Q + 0, it yields 

\d(xz ■ Ktfe- 2 "^ < n 2 \5xz(Q\ 2 e 2nQ{w) - cn6 ". 
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We have shown that 

P(xz--K?J <ne- nc6 «e nQ{w \ 

II V 'WnQ 

In view of the estimate (|4.11t , we then have 

I|m||„q < ^ie- nal e nQ(w) . 

Since u is analytic in D(z; 1/ yfn) we can now invoke the simple estimate (e.g. J2], Lemma 3.2) 

|fi(z)|V 2n( 2( z ) < n\\u\\ 2 nQ 

to get 

\u(z)\ < Me -»<^ e tt(Q(z)+Q>))_ 
This gives (|4.8b , and finishes the proof of Lemma l4~4l □ 
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